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We study strong tunneling (i.e. transmission h/e 2 Rt 3> 1) 
in the single-electron box with many transverse modes at 
zero temperature. We develop a new renormalization group 
method which includes all charge states and requires no ini- 
tial or final energy cutoff. We determine the ground-state 
energy, the average charge and the renormalized charging en- 
ergy. The covered range for the coupling constant and the 
gate voltage is much increased in comparison to recent per- 
turbative approaches, poor man scaling methods and Quan- 
tum Monte Carlo simulations. We reach the regime where 
Coulomb blockade become practically unobservable. 
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Introduction. Metallic nanostructures with junctions 
of low capacitance show quantization of electric charge 
and Coulomb blockade phenomena due to large charg- 
ing energy The simplest system in which one can 
see Coulomb blockade effects is the single-electron box. 
Electrons enter or leave a small metallic island via a tun- 
nel junction with many transverse channels. The junc- 
tion is characterized by a resistance Rt and a capaci- 
tance C . The charge on the island is controlled by a 
capacitively (with Cq) coupled external voltage Vq. In 
the "weak tunneling" regime, i.e., when the transmis- 
sion T = h/e 2 RT of the barrier is much less than unity, 
the average charge of the box can be described within 
the "orthodox theory" (!]] which treats tunneling in low- 
est order perturbation theory (golden rule). This means 
that the energy for charge N on the island is given by 
the classical expression 

E N (n x ) = E c (N-n x ) 2 , (1) 

with E c = e 2 /2(C + C G ) and n x = C G V G . At zero 
temperature, the ground-state energy as function of n x 
is the periodic continuation of the parabola Eo(n x ) for 
— 1/2 < n x < 1/2. Furthermore, the average charge, 
given by the periodic continuation of 

1 dE Q 



(3) 



-(nx) , 



(2) 



2E C dn x 

increases in steps at half-integer values of n x . 

In this letter, we discuss the fundamental question how 
quantum fluctuations influence Coulomb blockade phe- 
nomena at zero temperature, i.e., we analyze the smear- 
ing of E (n x ) and (N)(n x ) in the whole regime from weak 
to strong tunneling up to T ~ 50 (in recent experiments 
T ~ 33 was achieved Q]). Another quantity of recent 
interest is the second derivative of the energy Eq, 



Perturbative approaches in a — T/(47r 2 ) up to a 3 
H can describe E c (0) up to T ~ 15. However, within 
this approach, the range of allowed values for T decreases 
considerably for n x 1/2, and the average charge turns 
out to be divergent at n x = 1/2. 

Therefore, nonperturbative approaches have been de- 
veloped to describe the system near the degeneracy point 
of adjacent charge states either within a "non-crossing" 
expansion scheme Q or using an effective two-state sys- 
tem -Q. But for T > 1 the broadening h/(R T C) ~ 
TEc/k of the charge states is of the same order as the 
energy gap E Gl i.e., all charge states become important. 
And even for weaker coupling the results from a renor- 
malization group (RG) treatment (^Jl^l depend on both 
a high-energy and a low-energy cutoff. The first one, of 
the order of the charging energy, accounts for the ne- 
glected higher charge states. The second one is provided 
by the dominant low-energy scale, energy gap or temper- 
ature. As a result, quantitative statements depend on 
adjustable parameters. 

For very strong coupling, T> 1, instanton methods 
have been developed, which predict an exponential de- 
crease of Coulomb blockade phenomena: E c (0)/Ec ~ 
T v exp (— T/2), jl3],[l4]]. However, the exponent n is still 
controversial. Poor man scaling methods in 1/T reveal 
a strong renormalization of T to smaller values ]l0| , p^ |. 
This analysis is performed in the phase representation for 
the topological sector with winding number zero. How- 
ever, for decreasing T all winding numbers become im- 
portant. Consequently, a quantitative matching to the 
two-charge state approximation is not possible |l6| . 

Even recent quantum Monte Carlo (QMC) simulations 
§,|5[0@ have led t0 different results for E c (0) in the 
regime T > 15. Furthermore, the error bars and the 
computing time are very large for T>1 and/or n x — > 
1/2 and extrapolation to zero temperature in order to 
describe ground state properties becomes difficult. As a 
consequence, a large gap still remains between the small 
and large T-regions, the latter being even controversial. 

It is the purpose of this letter to present an indepen- 
dent, real-time RG-method which covers self-consistently 
level-renormalization and -broadening together with ver- 
tex corrections, and includes all charge states without the 
need of initial or final energy cutoffs. In contrast to QMC 
simulations our theory works also at zero temperature. 
For E c (0) we find good agreement with the QMC sim- 
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ulation of Refs. Pjlq]. We also find agreement with the 
exponential behavior of ££(0) for T > 20, as predicted 
by the controversial instanton results in Refs. jl3],[l4| but 
with a different preexponent. We calculate (N)(n x ) for 
all n x and predict that for T > 27 charge quantization 
will be suppressed such that the deviation from a straight 
line is less than 1%. We find that the slope of (N) at 
n x — 1/2 is always infinite. 

Model. The system is modeled by the Hamiltonian 
H = Hq + Ht , in which 

H o=Yl e kn a\ n a kn + £qnc\ n c qn + ^ E N \N)(N\ 



N 



(4) 



describes the box without tunneling. The operators a^) 
and c^J annihilate (create) electrons in the lead and is- 
land electrode, respectively. The tunneling coupling is 
given by 



H T = J2J29N\N + a){N\J a 



(5) 



N <r=± 



with J_ = J2kqn T k q a kn c gn and J + = JL . The channel 
number n = 1, . . . , N c h accounts for both the transverse 
wave number and the spin and is conserved during a tun- 
nel process. We assume here a wide metallic junction 
with a large number of transverse channels. The matrix 
elements TJ} q are assumed to be independent of the states 
k and q. They are related to the tunnel resistance by 
1/R T = (2^e 2 /fr)E„^(0)iV/(0)|T™| 2 , where N L/I (0) 
are the density of states of the lead and the island. The 
factor describes the relative strength of the tunnel 
coupling. At the beginning, we have g^ = 1 for all N 
and a. But under the renormalization group procedure, 
the coupling will indeed depend on the charge state N 
and on a. We define = ag^ f g^ !+1 . 

Renormalization group procedure. The invari- 
ant quantity under the RG is the S-matrix S = 
Texp(-i J dtH{t)). The idea of the RG is to start with 
some high-energy cutoff, reducing it by integrating out 
the corresponding degrees of freedom, and end up with an 
effective low-energy theory. All the energies and coupling 
constants will be renormalized, and, simultaneously, new 
terms which are not present in the original Hamiltonian 
may be generated. We will proceed in three steps, (i) 
We introduce a time cutoff t c . (ii) The cutoff is increased 
t c — ► t c + dt c , thus flowing from t® c to t[. The invariance 
condition for the S-matrix demands 



Te 



-i J dtH(t) I 



Te~ 



i J dtH'(t) I 



\t c +dt c 



(0) 



for each RG step, (iii) We take t Q c — > and t[ — > oo. 

(i) The S-matrix reads in interaction picture S = 
Texp(-i / dtH )T exp(-i fdtH T (t)i) in which H T (t)i 



is the tunnel part of the Hamiltonian in interaction pic- 
ture with respect to Hq. The expansion of the second 
exponential yields 



5' 



n=0 



>V I H x ...H n , (7) 

t 1 >t 2 >—>t n 



where Hi = Ht(U)i- Next, we perform a normal or- 
dering of the fermion operators c^) and involved in 

Hi, 

H l = :H 1 : (8) 

H X H 2 = : HiH 2 : +Hjl 2 (9) 
H\HiH-$ = : H1H2H3 : +H1H2H3 

+HxH 2 H s + fTifl^ffa (10) 

where : . . . : stands for normal ordering and H\H 2 repre- 
sents contractions of all the Fermi operators and 
in Hi and H 2 with the relative time t\ — t% ■ In principle 
the two Fermi operators of Hi may be contracted with 
Fermi operators of different Hi^s. For a large number 
of channels, N c ^ — > 00, however, these contributions are 
negligible. The only combination of contractions which 
enter here is (at zero temperature) 



(J±(h)J T (t 2 )) eq 



+ iair6'{tx-t 2 ). (11) 



We now introduce a cutoff t c by multiplying the right 
hand side of Eq. (JlT]) with 0(|ti — t 2 \ — t c ). Conse- 
quently, via Eqs. (@)- (|l(]) , the S-matrix is well-defined 
with a sharp cutoff t c . 

(ii) The increase of the cutoff t' c = t c + dt c leads to cor- 

1 — * — 1 

rcction terms H\H 2 , H\H 2 H-i, etc. The cross indicates 
that the time difference is between t c and t' c . We inter- 
pret such many-time objects as clusters with one time 
argument for the time ordering (we choose t 2 ). ft is pos- 
sible to reexponentiate the sum of terms and correction 
terms. The renormalized exponent reads 



Hi 



ffi + (-») 



rX-i 

H\H 2 



h>t 2 



( 1 — X — 1 r*i -i 

<^HiH 2 H 3 — H 2 HiH 3 j (12) 

tl>*2>*3 

plus higher-order correction terms which we truncate. 
This is the central equation since it determines the renor- 
malization of all quantities. The second term (propaga- 
tor renormalization) can be absorbed in a change of the 
energies, while the third term (vertex corrections) renor- 
malizes the coupling strength. 

(iii) The limit t° — > is possible if we take into account 
all charge states. There are no divergences like in a two- 
state approximation. After performing the limit t[ — > 00 
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we end up with an effective theory for the charge degrees 
of freedom without any contractions. 

To present the results, we introduce the abbreviations 
/ = ln(f c /t°), A N = E N+1 -E N , E N = E N t c , and A N = 
Ajv^c- The RG equations, then, are 



dE 



N 



dl 
dA N 
dl 



A N + i ( a N e lAN + aN+ie 



n -iA» 



da 



N 



,— iA» 



ctN-ie 



iAjv-l 



(13) 



(14) 



dl 



-i(Ajv+i - A N ) 



+2a N -i- 



giAjy _ giAjv-i 



i(A N - AjV_l) 

-a N e lKN - a N+1 e~ lAN + 1 



-a.Ne 



-iAj- 



a>N-ie 



„iAjv_ 



'} 



(15) 



Discussion of the RG equations. We note the follow- 
ing: 

- All energies become complex, which is related to the 
broadening of the levels. 

- The coupling strength apj and the difference A^v+i — 
Ajv, which are initially a and 2Eq , respectively, evolve 
different for different charge states N. 

- The RG equations contain all terms from perturbation 
theory up to Q(a 2 ). 

- Eqs. ( |l4|) and ( |l5|) include also all terms of the poor 
man's scaling. If one considers for n x m 1/2 only two 
charge states N — 0,1, and expands all exponentials 
exp(±iA ) Ml± iA Q , Eqs. (0) and @ reduce to the 
well-known formulas dAo/dl = Ao(l— 2ao) and da^/dl = 
-2a\, with solutions A (Z)/A (0) = a (0/ a o(0) = 
[1 + 2a (0)Zp 1 with a (0) = a. 

Results. In Figs. |l|||, we show Bq, and Eq as 

function of n x up to a = 0.5 (T ~ 20) and compare with 
QMC data from Refs. P,fl8||. (In the numerics we had 



to choose the initial cutoff much smaller than O.IE^ 1 . 
Here, we chose O.OlEg 1 .) As expected we find a grad- 
ual flattening of the ground-state energy with increasing 
coupling and a rounding of the average charge towards 
a linear function. The average charge (N)(n x ) agrees 
well with available QMC data for a < 0.3, except for 
n x ^ 0.4. This deviation can be traced back to a finite 
temperature used in the QMC simulations ^From 
our results we conclude: (i) For T ~ 27, the shape of 
(N)(n x ) is nearly indistinguishable from a straight line 
within an error of 1%. (ii) The slope of (N) at n x — 1/2 
is infinite, i.e., there is always a very small crossover re- 
gion around the degeneracy point showing a deviation 
from a straight line. This is explicitly shown in Fig. ^ 
for the renormalized charging energy which diverges at 
n x = 1/2. However, this region is so small that it can 



hardly be identified experimentally, (iii) The renormal- 
ized charging energy at n x = decreases with increasing 
coupling. This reveals another indication for the washing 
out of charge quantization and is shown in more detail 
in Fig. ^ using a logarithmic scale for Eq(0) as function 
of the coupling a. Here we find very good agreement of 
our results with perturbation theory |^],[3 up to a ~ 0.4 
(T - 15) and with the QMC data of Ref. [§. Fig. | shows 
the comparism to other QMC simulations. Including er- 
ror bars, our data are consistent with the QMC data of 
Ref. jD| (the authors report that the saturation for the 
largest T is still unclear [^0| ) , but are slightly lower than 
those of Ref. , and higher than those of Ref. @ . For 
T > 25 no QMC data are available yet but the numeri- 
cal solution of the RG-equations (|l3|)-([l5|) is stable until 
T ~ 50. For T > 20 our data are consistent with the 
exponential behavior E^(0)/E c ~ T ?7 exp(-T/2) pre- 
dicted by instanton techniques but with i] = 6.5 
rather than n — 2 or 3. 

Summary. We have investigated strong tunneling in 
the single-electron box with the help of a new, renormal- 
ization group method beyond the perturbative regime. 
All charge states are included and, therefore, all results 
are cut-off independent. We analyzed the observability 
of charging effects at zero temperature by computing the 
average charge as function of the gate voltage. We found 
a deviation of less than 1% from a straight line for T > 27 
and an infinite slope at the degeneracy point. Further- 
more we compared our results to contradicting QMC sim- 
ulations and instanton techniques. 
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FIG. 1. Ground-state energy as function of n x . 
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FIG. 2. Average charge as function of n x . The step at 
n x — 1/2 is rounded for increasing a, approaching more and 
more a straight line. The circles are QMC data from Ref. |l£ 
for a = 0.3, 0.2, 0.1. The triangles are QMC data from Ref. 
for a = 0.25330, 0.12665. 
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FIG. 3. Eq as a function of n x . For n x — > 1/2 it diverges. 
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FIG. 4. Eq as function of a for n x = 0. Solid 
line: renormalization group method described in this let- 
ter. Dashed lines: perturbation theory of order a, a 2 , 
and a 3 from Refs. lH,q|. Long-dashed lines: large-a expan- 
sions from Refs. |l4]| (upper curve) and ^] (lower curve). 
Dot-dashed line: For a > 0.5 our data are consistent with 
0.5(27ra) 6J5 exp(-27r 2 a). Circles: QMC data from Ref. §]. 
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FIG. 5. Eq as function of a for n x = 0. Solid line: renor- 
malization group method described in this letter. Dashed 
lines: perturbation theory of order a 3 ||. QMC data from 
Ref. @ (o), § (o), H '(□) and g (*). 
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